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Cancellation of laser frequency noise in interferometers is crucial for attaining the requisite 
sensitivity of the triangular 3-spacecraft LISA configuration. Raw laser noise is several orders of 
magnitude above the other noises and thus it is essential to bring it down to the level of other 
noises such as shot, acceleration, etc. Since it is impossible to maintain equal distances between 
spacecrafts, laser noise cancellation must be achieved by appropriately combining the six beams 
with appropriate time-delays. It has been shown in several recent papers that such combinations 
are possible. In this paper, we present a rigorous and systematic formalism based on algebraic 
geometrical methods involving computational commutative algebra, which generates in principle all 
the data combinations cancelling the laser frequency noise. The relevant data combinations form 
the first module of syzygies, as it is called in the literature of algebraic geometry. The module is 
over a polynomial ring in three variables, the three variables corresponding to the three time-delays 
around the LISA triangle. Specifically, we list several sets of generators for the module whose linear 
combinations with polynomial coefficients generate the entire module. We find that this formalism 
can also be extended in a straight forward way to cancel Doppler shifts due to optical bench motions. 

The two modules are infact isomorphic. 

We use our formalism to obtain the transfer functions for the six beams and for the generators. 

We specifically investigate monochromatic gravitational wave sources in the LISA band and carry 
out the maximisiation over linear combinations of the generators of the signal-to-noise ratios with 
the frequency and source direction angles as parameters. 

I. INTRODUCTION 

Breakthroughs in modern technology have made possible the construction of extremely large interferometers both 
on ground and in space for the detection and observation of gravitational waves (GW). Several ground based detectors 
are being constructed around the globe; these are the projects, LIGO, VIRGO, GEO, TAMA and AIGO of building 
interferometers whose armlengths will be of the order of kilometers. These detectors will operate in the high frequency 
range of GW of ~ 10 Hz to a few kHz. A natural limit occurs on decreasing the lower frequency cut-off of 10 Hz 
because it is not practical to increase the armlengths on ground and also because of the gravity gradient noise which 
is difficult to eliminate below 10 Hz. Thus, the ground based interferometers will not be sensitive below the limiting 
frequency of ~ 10 Hz. But on the otherhand, there exist in the cosmos, interesting astrophysical GW sources which 
emit GW below this frequency such as the galactic binaries, massive and supermassive blackhole binaries etc. If 
we wish to observe these sources, we need to go to lower frequencies. The solution is to build an interferometer in 
space, where such noises will be absent and allow the detection of GW in the low frequency regime. LISA - Laser 
Interferometric Space Antenna - is a proposed mission which will use coherent laser beams exchanged between three 
identical spacecrafts forming a giant (almost) equilateral triangle of side 5 x 10 6 kilometers to observe and detect low 
frequency cosmic GW. The ground based detectors and LISA complement each other in the observation of GW in 
an essential way, analogous to the optical, radio, X-ray, y-ray etc., observations do for the electromagnetic waves. As 
these detectors begin to operate, a new era of Gravitational Astronomy is on the horizon and a radically different 
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view of the universe is expected to be revealed. 

In ground based detectors the arms are chosen to be of equal length so that the laser light experiences identical 
delay in each arm of the interferometer. This arrangement precisely cancels the laser frequency/phase noise at the 
photodetector. This cancellation of noise is crucial since the raw laser noise is orders of magnitude larger than 
other noises in the interferometer. The required sensitivity of the instrument can thus only be achieved by near 
exact cancellation of the laser frequency noise. However, in LISA it is impossible to achieve equal distances between 
spacecrafts and the laser noise cannot be cancelled in an obvious manner. In LISA, six data streams arise from the 
exchange of laser beams between the three spacecrafts - it is not possible to bounce laser beams between different 
spacecrafts, as is done in ground based detectors; a scheme analogous to the RF transponder scheme is used as was 
done in the early experiments for detecting GW by Doppler tracking a spacecraft from Earth. Several schemes, some 
quite elaborate, have been proposed [|],[|] which combine the recorded data by suitable time-delays corresponding 
to the three armlengths of the giant triangular interferometer. These schemes have data combinations of six or at 
most eight data points which give respectively a six and eight pulse response of GW and also show how other data 
combinations can be obtained by linear superposition. 

Galactic and extragalactic binaries are important sources in the LISA frequency band. Their abundance and 
resulting spectral amplitude has been estimated using the population synthesis by various authors ill- In the lower 
frequency range (< 1 mHz) there are a large number of such sources in each of the frequency bins. This makes it 
impossible to resolve an individual source which results in a stochastic GW background. It has been also proposed that 
the galactic halo MACHOs such as white dwarfs and blackholes ( with mass ~ 0.5 M 0 ) can also produce stochastic 
GW background pli[|T8| . In a recent work, Tinto at. al |lij have used Doppler delayed beams for discriminating the 
stochastic background from the instrumental noise. The angular resolution of the LISA is restricted because it is an 
all-sky monitoring detector with quadrupole beam pattern, however, the angular resolution can be achieved by the 
relative amplitude and phases of the two polarisations and Doppler modulation of the beams due to the motion of 
LISA around the sun j2(j,[2jj]. 

We start with the fundamental papers by the JPL team 01 where idea of delayed data combinations was first 
proposed. Here we present a systematic method based on modules over polynomial rings, which not only reproduces 
the previously obtained results, but guarantees all the data combinations which cancel the laser noise. The data 
combinations in the case of laser frequency noise consist of the six suitably delayed data streams, the delays being 
integer multiples of the light travel times between spacecrafts, which can be conveniently expressed in terms of 
polynomials in the three delay operators Ei, E 2 , E$ corresponding to the light travel time along the three arms. The 
laser noise cancellation condition puts three constraints on the six polynomials of the delay operators corresponding to 
the six data streams. The problem therefore consists of finding six tuples of polynomials which satisfy the laser noise 
cancellation constraints. These polynomial tuples form a modulc0 called in the literature, the module of syzygies. 
There exist standard methods for obtaining the module, by that we mean, methods for obtaining the generators of 
the module so that the linear combinations of the generators generate the entire module. Three constraints on six 
tuples of polynomials do not lead to three generators as naive reasoning might suggest. Here we are dealing with 
modules rather than vector spaces and the rules are different. The procedure first consists of obtaining a Groebner 
basis for the ideal generated by the coefficients appearing in the constraints. This ideal is in the polynomial ring 
in E\, E 21 E$ over the domain of rational numbers (or integers if one gets rid of the denominators). To obtain the 


1 A module is an abelian group over a ring as contrasted with a vector space which is an abelian group over a field. The scalars 
form a ring and just like in a vector space, scalar multiplication is defined. However, in a ring the multiplicative inverses do 
not exist in general for the elements, which makes all the difference! 
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Groebner basis for the ideal, one may use the Buchberger algorithm or use a package such as Mathematica. From 
the Groebner basis there is a standard way to obtain a generating set for the required module. All of this procedure 
has been described in the literature [^,^|. We thus obtain seven generators for the module. However, the method 
does not guarantee a minimal set and we find that a generating set of 4 polynomial six tuples suffice to generate 
the required module. Alternatively, we can obtain generating sets by using the software Macaulay 2. It gives us a 
a Groebner basis for the module consisting of five generators and another generating set consisting of six elements. 
The importance of obtaining more data combinations is evident: they provide the necessary redundancy - different 
data combinations produce different transfer functions for GW and so specific data combinations could be optimal 
for given astrophysical source parameters in the context of maximising signal-to-noise (SNR), detection probability, 
improving parameter estimates etc. 

The scheme we have described above can also be extended in a straight forward way to include optical bench 
motions. There are now twelve Doppler streams of data and we apply the above scheme to cancel the noise due to 
optical bench drift and laser frequency noise. The six extra streams can be combined two by two by subtracting one 
stream from the other to obtain three streams in which the frequency shifts in the optical fibers are cancelled. Thus 
we have only nine streams to contend with and now the module consists of nine tuples of polynomials on which six 
linear constraints are imposed. We show that the problem can be solved in the terms of the previous one where the 
three extra polynomials are written in terms of the six tuple polynomials which are solutions to the laser frequency 
noise cancellation problem. Thus the solution to the first problem extends easily to the second. 

Finally, we apply our formalism to a class of important astrophysical sources, but relatively simple to analyse, 
namely, monochromatic GW sources. We maximise the SNR for such sources over much of the module of data 
combinations by considering linear combinations of the generators of the module with the coefficients being real 
numbers. Strictly speaking one must take polynomials as the coefficients so that the maximisation extends to the 
entire module, but we find that even this simplifying assumption yields satisfactory results. We present the maximised 
SNR as a function of frequency over the data combinations. 

We organise the paper in the following manner: In the section II, we present the six raw data streams obtained with 
the laser phase noise and formulate the conditions for the laser phase noise cancellation. We also obtain difference 
equations which should be satisfied by the time-delay operators for cancelling the laser noise. The solutions for the 
noise cancellation combinations can be represented as the modules syzygies over the polynomial ring using standard 
methods of algebraic geometry described in the section III. First a Groebner basis for the ideal is obtained. From 
the Groebner basis the generators for the module of syzygies can be computed. Several sets of generators have been 
obtained for this module. In the subsection III B this approach is extended to cancel the acceleration noise of the 
optical benches. In the section IV we compute the detector response for the GW signal and obtain transfer functions 
for the six elementary beams. In the section V, first, we determine the effective noise for each generator by taking 
shotnoise and acceleration noise of the proof masses into account, which do not cancel out in the combinations. We 
obtain SNRs for monochromatic sources and maximise the SNR over the allowed data combinations that cancel the 
laser frequency noise. 


II. TIME DELAYED DATA AND THE DIFFERENCE EQUATION 

A. Time-delayed Data 

We label the spacecrafts as 1,2 and 3. Let L\, L 2 , L 3 be the lengths of the arms (sides of the triangle) where L 3 is 
the distance between spacecrafts 1 and 2; and so on by cyclic rotation of indices (see figure |l]). Each spacecraft has a 
laser which is 
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• used to send beams to the other two spacecrafts, and, 

• used as a local oscillator to produce a beat signal with the incoming beams from the other two spacecrafts. 


The data are recorded as fractional Doppler shifts. These fractional Doppler shifts can occur due to the GW signal 
and the noise. Here we will be concerned with the laser frequency noise only. More precisely, if vq is the central 
frequency of the laser and the frequency fluctuation of the laser on spacecraft i at time t is Ag; (t ), then the fractional 
frequency fluctuation Ci(t) is given by, 


aft) 


A Uj(t) 
v 3 


( 2 . 1 ) 


The six streams of Doppler data are obtained from the Ci(t) by combining them suitably with their time delayed 
copies, where the time delays are just the light travel times between the three spacecrafts. We adopt the following 
notation for the six streams: we divide the six streams into two sets U l and V 1 where, i = 1,2,3 of three streams 
each. U l and V 1 can be regarded as 3 component vectors U and V respectively. Also we denote the time-delayed 
data in arm k, k = 1,2,3, by the shift operator E k whose action on a function f(t) is described by the equation: 



FIG. 1. LISA Geometry 

E k f(t)= f(t-L k ), (2.2) 

where we have chosen units in which the speed of light is unity. The six streams are: 

u 1 = E 2 C 3 - Cl, 

U 2 = E 3 Ci — C2 , 

U 3 = E\C 2 — C.3, 

V 1 = C\ — E 3 C 2 , 

V 2 = C 2 - E 1 C 3 , 

V 3 = C 3 -E 2 a. (2.3) 

Thus explicitly we have, t/ 1 (t) = C 3 (t — L 2 ) — C\{t) etc. t/ 1 (t) is the data stream obtained by beating the laser beam 
transmitted by spacecraft 3 to spacecraft 1 measured at time t at spacecraft 1; and so on by cyclic rotation. Similarly 
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— V 1 ^) = C 2 (f — L 3 ) — Ci(t) is the laser beam emitted by spacecraft 2 and received and beaten with the laser at 
spacecraft 1 at time t. If we label the spacecrafts in a counter-clockwise (positive sense) fashion, then the beams U l 
travel in the positive sense while the beams V x travel in the negative sense. 

The goal of the analysis is to add suitably delayed beams together so that the net result is zero. This amounts to 
seeking data combinations which cancel the laser frequency noise. In the notation/formalism that we have invoked, 
the delay is obtained by applying the operators Ek to the beams U l and V 1 . A delay of k\Li + fc 2 L 2 + k^L^ is 
represented by the operator E^E^Eff acting on the data, where k\ , fc 2 and are integers. In general a polynomial 
in Ek, which is a polynomial in three variables, applied to say U 1 combines the same data stream C7 1 (i) with different 
time-delays of the form k\L\ + fc 2 L 2 + /C 3 L 3 . This notation conveniently rephrases the problem. One must find six 
polynomials say pt{E\, E 2 , £ 3 ), < 7 j(£i, £ 2 , £ 3 ), i = 1, 2,3 such that: 

3 

Y,PiV i +q i U i =0. (2.4) 

i=i 

Cancellation of the laser noise is implicit in the above equation. 


B. The Difference Equation for Shift Operators 


It is convenient to express eq. (2.3) in matrix form. This allows us to obtain a matrix operator equation whose 
solutions are p and q where we have now written p l and q l as column vectors. We can similarly express U l , V 1 , C l as 


column vectors U, V, C respectively. In matrix form eq. (2.3) become: 


V = /r T C, 


U = -n C, 


where, /1 is a 3 x 3 matrix given by, 


T = 



The exponent ‘T’ represents the transpose of the matrix. The eq. (2H) becomes: 

(mp - M T q) T C = 0 . 


(2.5) 


( 2 . 6 ) 


(2.7) 


where we have taken care to put C on the right of the operators. Since the above equation must be satisfied for 
arbitrary ‘data’ C, we obtain a matrix equation for the shift operators: 


mp = fi T q- 


( 2 . 8 ) 


The solutions to eq. (2.8) are p, q which are column vectors of polynomials in the shift operators Ek- Note that since 
the Ek are just shift operators, they commute, so the order in writing these operators is unimportant. In mathematical 
terms, the polynomials form a commutative ring. 

We can formally solve for p since the matrix /r is invertible. However, det p, = I—E 1 E 2 E 3 appears in the denominator 
on the R.H.S., which leads to the division by polynomials in Ek- This may seem hard to interpret. But we can pull 
this factor to the L. H. S. to ‘rationalise’ the expressions. Then we obtain, 


Ap = A q, 


(2.9) 


where A = n a djP T and p. a dj is the adjoint of /i. The operator A = 1 — E\ E 2 E 3 is the usual difference operator that 
appears in finite differences and difference equations. The quantity E\ E 2 E 3 plays a central role in determining the 
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natural time-step for the problem, namely, s = Li + L 2 + £ 3 ; which is nothing but the light travel time around the 
perimeter of the LISA triangle. A is just the difference corresponding to this time-step. 


Explicitly, using (2.6) the matrix A is given by: 


A = 


/ 1 -El E\E'2 - E :i E. 2 (l - El) 

(1 - E'i)E :i 1 - El E 2 E 3 - E 1 
\ E 3 Ei — E 2 Ei{\ - E 3 ) l-El 


( 2 . 10 ) 


The equations display a cyclic symmetry in the indices 1,2,3 which is also apparent in the matrix A. The cyclic 
symmetry results from the nature of the problem since we are free to choose the labeling of the three space-crafts. In 
the matrix A we must also change the rows/columns consistently performing the cyclic change of the indices. The 
cyclic symmetry is further carried over to the solutions (p, q). 


The integration of the eq. (T9) can be carried out in time-steps of s. The integration is immediate if we operate 


the eq. (T9) on V. We first need to take the transpose of eq. (2J) and then operate on V. We then obtain: 

Ap T V = q T A T V 

= q T (/WiM T ) T v 
= q T fj, A C 

= -Aq T U, (2.11) 


which gives, 


A(p T V + q T U) = 0. 


( 2 . 12 ) 


This equation immediately integrates to, 

(p T V + q T U)(f + ns) = (p T V + q T U)(f), (2-13) 

where n is an integer. If we arbitrarily set t = 0 and if (p T V + q 7 U)(0) = 0, then (p T V + q T U)(ns) = 0. 

It is not clear to us, how the above solution would be useful physically, but we present it as an interesting outcome. 


However, the main problem is of seeking solutions (p, q) to eq.(2.8). This problem and its solution we discuss in the 
next section. 


III. THE MODULES OF SYZYGIES 


Several solutions have been exhibited in |l], |2j to eq.(T 8 ). The solutions have the characteristic property that 
the det fi cancels on both sides leading to polynomial vectors p and q. We reproduce here the solutions obtained in 
previous works. The solution £ is given by p T = q T = (E\,E 2 ,E 3 ). The solution a is described by p T = (1,E 3 ,EiE 3 ) 
and q T = (1, EiE 2 , E 2 ). The solutions (3 and 7 are obtained from a by cyclically permuting the indices of Ek, p and 
q. These solutions as realised in earlier works are important, because they consist of polynomials with lowest possible 
degrees and thus are simple. Other solutions containing higher degree polynomials can be generated conveniently 
from these solutions. Linear combinations of these solutions are also solutions to the given system of equations. But 
it is not clear that this procedure generates all the solutions. In particular, £ cannot be generated from the set a, (3 
and 7 , where generating a data combination means writing it as a linear combination of the elements belonging to 
the generating set. 

The basic reason, as mentioned earlier, is that we do not have a vector space here. Three independent constraints 
on a six tuple do not produce a space which is necessarily generated by three basis elements. This conclusion would 
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follow if the solutions formed a vector space but they do not. The polynomial six-tuple p, q can be multiplied by 
polynomials in £ 1 , £ 2 , £3 (scalars) which do not form a field. So that the inverse in general does not exist within the 
ring of polynomials. We therefore have a module over the ring of polynomials in the three variables E±, E 2 , £ 3 . 

In this section we present the general methodology for obtaining the solutions to (2.8). The method is illustrated by 
applying it to equations ( |2.£j ). In the next subsection we consider the more general problem of optical bench motions 
as well. The optical bench motion noise can also be eliminated using the same method. 


A. Cancellation of laser frequency noise 


There are three linear constraints on the polynomials given by the equations (T 8 ). Since the equations are linear 
the solutions space is a submodule of the module of six-tuples of polynomials. The module of six-tuples is a free 
module, i.e. it has six basis elements that not only generate the module but are linearly independent. A natural 
choice of the basis is fi = (0,..., 1,..., 0) with 1 in the i-th place and 0 everywhere else; i runs from 1 to 6 . The 
definitions of generation (spanning) and linear independence are the same as that for vector spaces. A free module is 
essentially like a vector space. But our interest lies in its submodule which need not be free and need not have just 
three generators as it would seem if we were dealing with vector spaces. 

The problem at hand is of finding the generators of this submodule i.e. any element of the module should be 
expressible as a linear combination of the generating set. In this way the generators are capable of spanning the full 
module or generating the module. We examine the eq.(2.8) explicitly componentwise: 


Pi~qi + E 3 q 2 - E 2 p 3 = 0 , 

P2 - 92 + £i93 ~ E 3 pi = 0, 

P 3 ?3 t £29l E1P2 — 0. (3.1) 


The first step is to use Gaussian elimination to obtain pi and P2 in terms of p 3 , qi, 92, 93- We then obtain: 


Pi = qi~ E 3 q 2 + E 2 p 3 , 

P2 = 92 - £193 + E 3 pi 

= £391 + (1 — £3)172 ^ £i <73 + E2E3P3, ( 3 - 2 ) 

and then substitute these values in the third equation to obtain a linear implicit relation between 723,91,92,93- We 
then have: 


(■ E 1 E 2 E 3 — V)P 3 + {E 1 E 3 — £2)91 + £i(l — £3)92 + (1 — £1)93 = 0 . 


( 3 . 3 ) 


Obtaining solutions to (3.3) amounts to solving the problem since the the remaining polynomials Pi,P 2 have been 
expressed in terms of 723 , 91,92,93 in (3.2). 

The solutions to (3.3) form the first module of syzygies of the coefficients in the eq. ( |3.3[ ), namely, 

E 1 E 2 E 3 — 1, E 1 E 3 — £ 2 , £i(l — £|), 1 — Ef. The generators of this module can be obtained from standard methods 
available in the literature. We briefly outline the method given in the books by Becker et al. 0 and Kreuzer and 
Robbiano [p| below. The details have been included in the Appendix. 


1. Groebner basis 

The first step is to obtain the Groebner basis for the ideal U generated by the coefficients: 
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U\ — E1E2E3 — 1 , U2 — E1E3 — E2, U3 — Ei(l — E 2 ), 


U4 = 1 — E\. 


(3.4) 


The ideal IA consists of linear combinations of the form y) VjUj where Vi, i = 1,4 are polynomials in the ring 
Q[Ei, E 2 , E 3 ] where Q is the field of rational numbers. There can be several sets of generators for U. A Groebner 
basis is a set of generators which is ‘small’ in a specific sense. 

There are several ways to look at the theory of Groebner basis. One way is, suppose we are given polynomials 
51, <72, • 5m in one variable over say Q and we would like to know whether another polynomial / belongs to the ideal 
generated by the g’s. A good way to decide the issue would be to first compute the gcd (greatest common divisor) 
g of 51,52, 5m and checking whether / is a multiple of 5. One can achieve this by doing the long division of / 
by 5 and checking whether the remainder is zero. All this is possible because Q[x] is a Euclidean domain and also 
a principle ideal domain (PID) wherein any ideal is generated by a single element. Therefore we have essentially 
just one polynomial - the gcd - which generates the ideal generated by 51,52, ...,5m- The ring of integers or the ring 
of polynomials in one variable over any field are examples of PIDs whose ideals are generated by single elements. 
However, when we consider more general rings (not PIDs) like the one we are dealing with here, we do not have a 
single gcd but a set of several polynomials which generates an ideal in general. A Groebner basis of an ideal can be 
thought of as a generalisation of the gcd. In the univariate case, the Groebner basis reduces to the gcd. 

Groebner basis theory generalises these ideas to multivariate polynomials which are neither Euclidean rings nor 
PIDs. Since there is in general not a single generator for an ideal, Groebner basis theory comes up with the idea of 
dividing a polynomial with a set of polynomials, the set of generators of the ideal, so that by successive divisions by 
the polynomials in this generating set of the given polynomial, the remainder becomes zero. Clearly, every generating 
set of polynomials need not possess this property. Those special generating sets that do possess this property (and 
they exist!) are called Groebner bases. In order for a division to be carried out in a sensible manner, an order must 
be put on the ring of polynomials, so that the final remainder after every division is strictly smaller than each of 
the divisors in the generating set. A natural order exists on the ring of integers or on the polynomial ring Q(x)', 
the degree of the polynomial decides the order in Q(x). However, even for polynomials in two variables there is no 
natural order apriori ( Is x 2 + y greater or smaller than x + y 2 ? ). But one can, by hand as it were, put a order on 
such a ring by saying x » y, where >> is an order, called the lexicographical order. We follow this type of order, 
Ei » E 2 » E 3 and ordering polynomials by considering their highest degree terms. It is possible to put different 
orderings on a given ring which then produces different Groebner bases. Clearly, a Groebner basis must have ‘small’ 
elements so that division is possible and every element of the ideal when divided by the Groebner basis elements 
leaves zero remainder, i.e. every element modulo the Groebner basis reduces to zero. 

In the literature, there exists a well-known algorithm called the the Buchberger algorithm which may be used to 
obtain the Groebner basis for a given set of polynomials in the ring. So a Groebner basis of IA can be obtained 


from the generators u t given in eq.(3.4) using this algorithm. It is essentially again a generalisation of the usual long 
division that we perform on univariate polynomials. More conveniently, we prefer to use the wellknown ‘Mathematica’ 
package. Mathematica yields a 3 element Groebner basis Q for U\ 


g = {E 2 -l,El-l,Ei~E 2 E 3 }. 


(3.5) 


One can easily check that all the it, of eq.(T4) are linear combinations of the polynomials in Q and hence Q generates U. 
One also observes that the elements look ‘small’ in the order mentioned above. However, one can satisfy oneself that Q 
is a Groebner basis by using the standard methods available in the literature. One method consists of computing the 
S-polynomials (see Appendix A) for all the pairs of the Groebner basis elements and checking whether these reduce 
to zero modulo Q. 

This Groebner basis of the ideal U is then used to obtain the generators for the module of syzygies. 
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2. Generating set for the module of syzygies 


The generating set for the module is obtained by further following the procedure in the literature HI- The de¬ 
tails are given in the Appendix specifically for our case. We obtain 7 generators for the module. These generators 
do not form a minimal set and there are relations between them; in fact this method does not guarantee a mini¬ 
mum set of generators. These generators can be expressed as linear combinations of cc, /?, 7 , £ and also in terms of 
iW,!®,!! 3 ),!* 4 ) given below in eq.( |3.6| ). The importance in obtaining the 7 generators is that the standard 
theorems guarantee that these 7 generators do infact generate the required module. Therefore from this proven set of 
generators we can check whether a particular set is infact a generating set. We present several generating sets below: 

Alternatively, we may use a software package called ‘Macaulay 2 ’ which calculates the generators given the the 
equations (3T). Using Macaulay 2, we obtain six generators. Again, Macaulay’s algorithm does not yield a minimal 
set; we can express the last two generators in terms of the first four. Below we list this smaller set of four generators 
in the order A = (pi,p 2 ,p 3 , qi, <?2, <73): 


X (1) = (E 1 E 3 - E 2 , 0, El - 1,0, E 2 E 3 - E u El - 1), 

-X (2) = (E 1 ,E 2 ,E 3 ,E 1 ,E 2 ,E 3 ), 

X (3) = (1, -U 3 , E\E 3 , 1, E 3 E 2i -£/ 2 ), 

X (4) = {E 3 E 2 , 1 , E u E 3 , 1 , E 2 E 3 ). (3.6) 


Note that the last three generators are just X^ = a , A^ 4 ) = (3. But there is an extra generator X^ needed 

to generate all the solutions. 

Another set of generators which could be useful for further work is a Groebner basis of a module. The concept of a 
Groebner basis of an ideal can be extended to that of a Groebner basis of a submodule of {K[xi,x 2 ,..., x n ]) m where 
AT is a field, since a module over the polynomial ring can be considered as generalisation of an ideal in a polynomial 
ring. Just as in the case of an ideal, a Groebner basis for a module is a generating set with special properties. For 
the module under consideration we obtain a Groebner basis using Macaulay 2 : 

G (1) = (Ei,E 2 , E 3 , E 1 ,E 2 ,E 3 ), 

G {2) = (E x E 3 - E 2 , 0 , El - 1 , 0 , E 2 E 3 -E^El-l), 

G (3) = (E 1 E 2 ,1,E 1 ,E 3 ,1,E 2 E 3 ), 

G (4) = (1,E 3 ,E 1 E 3 ,1,E 1 E 2 ,E 2 ), 

G (5) = (E 3 (Ef - 1), 1 - El 0,0,1 - El E\(El - 1)). (3.7) 

Note that in this Groebner basis = ( = X^ 2 \ G^ = AlW, G^ = (3 = X^ 4 \ G^ = a = X^ 3 \ Only G^ is the 
new generator. 

Another set of generators are just a, / 3 , 7 and (. This can be checked using Macaulay 2 or one can relate a , f3, 7 
and £ to the generators A ( A \A = 1,2, 3,4 by polynomial matrices. In the Appendix, we express the 7 generators 
we obtained following the literature, in terms of a,/3 ,7 and (. Also we express a,/3, 7 and ( in terms of X^ A \ This 
proves that all these sets generate the required module of syzygies. 

The question now arises as to which set of generators we should choose which facilitates further analysis. The 
analysis is simplified if we choose a smaller number of generators. Also we would prefer low degree polynomials to 
appear in the generators so as to avoid cancellation of leading terms in the polynomials. By these two criteria we 
may choose, A ^ or a, (3, 7 , Among these two sets of generators, we arbitrarily make the choice of X^ A \ 
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B. Cancellation of noise from moving optical benches 


The work done in ]l],[|] can be conveniently reexpressed in our formulation and leads to further insights into the 
problem. 

There are two optical benches on each space-craft which have random velocities and are connected by optical fibers. 
The random velocities of the optical benches and the delay in the optical fibers are measured as further Doppler 
shifts apart from other noise and the GW signal. Since we are interested in the cancellation of laser frequency noise 
and motion of the optical benches, we write expressions for the beams containing only these quantities. The Doppler 
beams will of course contain other effects arising from shot noise, GW signal, motion of proof masses etc., but we will 
not write them in the expressions for the Doppler data because they are not relevant to the problem we are interested 
in. We follow the notations of The quantities pertaining to the left bench will be unstarred while that for the 

right bench are starred. There are now twelve Doppler data streams which need to be combined in an appropriate 
manner in order to cancel the noise from the laser as well as from the motion of the optical benches. The fractional 
frequency fluctuations of laser on the left optical bench i are denoted by and on the right optical bench i* by C *, 
the random velocities of the benches V^, V* and ft the frequency shifts in the optical fibers connecting the optical 
benches in space-craft i. We then have the following expressions for the four data streams pertaining to space-craft 1: 

u 1 = E 2 (C 3 - n 2 • V 3 ) - (n 2 • V| + Cl), 

V 1 = ~E 3 {C * 2 + n 3 • v*) + (Cr - n 3 ■ V,), 

Zi = C!~ C? +m- 2 n 3 • V!, 

z\ = C{ - Ci + 771 + 2n 2 • VI . (3.8) 

The other eight data streams on space crafts 2 and 3 are obtained by cyclic permutations of the indices in the above 
equations. Here n 2 denotes a unit vector in the direction from space-craft 1 to space-craft 3 and the remaining unit 
vectors n 3 and ni are obtained by cyclically permutating the indices. 

We find that the twelve Doppler data streams depend only on the particular combinations C\ — n 3 • Vi and 
Cl + n 2 • V); and their cyclic permutations. We define these combinations as Ci and Cl respectively, i.e.: 

C\ = C\ - n 3 ■ V | , 

Cl = Cl + n 2 ■ VJ , (3.9) 


and also their cyclic permutations. Then the expressions for the data streams simplify considerably. We write the 
expressions for the data streams corresponding to space-craft 1. Others are obtained as before by cyclic permutations. 

U 1 = E 2 C 3 - Cl , 

V 1 = ~E 3 C* 2 + Ch , 

Z l = \{z x -zl) 

= Ci - Cl. (3.10) 


The new variables Z l have been defined which automatically cancel the 7 ^. Also we 
same form as in eq. (2.3), except that the Ci are replaced by the Ci which account 
benches. 

The noise cancellation condition now becomes: 


note that the U 1 , V 1 have the 
for the motions of the optical 


P,r + % u l + nz i = o, 


(3.11) 
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where r, are polynomials in E\, E 2 , E 3 . Since the above equations should hold for any Ci, C*, we obtain the following 
equations the polynomials Pi,qi , r* must satisfy: 


pi + E 3 q 2 + n = 0, 

E 2 p 3 + qi + n = 0 , 

P2 + E 3 q 3 + r 2 = 0 , 

E3P1 +q 2 + r 2 = 0 , 

P3 + E 2 qi +r 3 = 0 , 

E 1 p 2 + q 3 + r 3 = 0 . (3-12) 


We now have a nine component polynomial vector. The solutions to ( p. 12| ) form another module of syzygies which 
is related in a simple way to the module obtained in just laser noise cancellation. Eliminating n from eq. ( 3.12| ) we 
obtain the same equations for pi and qi as in (|3T|). This enables us to extend the previously obtained generating sets 
to this module. Thus, thanks to the mapping of Ci, (C*) —* Ci, ( C'*), the two modules are isomorphic. We just state 
the remaining three entries (ri,r 2 ,r 3 ) of the generating sets keeping the same notation. The first set of 4 generators 
in the order (ri, r 2 , r 3 ) are: 


* (1) = (£ 2(1 - - El), 1 - El), 

X^ = —(E\ -f- E 2 E 3 , E 2 + E\E 3 , E 3 + E 3 E 2 ), 

X ^ = —(1 + E\E 2 E 3 , E\E 2 + E 3 , E\E 3 + E 2 ), 

I (4) =-(E 1 E 2 + E 3 ,l + E 1 E 2 E 3 ,E 1 +E 2 E 3 ). (3.13) 

In the other generating set, namely, the Groebner basis we need to specify just C 3 since the other elements are in the 
previous generating set. Thus, 


= ( 0 , (El - 1)(1 - El), 0 ). 


(3.14) 


IV. THE DETECTOR RESPONSE 

The ripples produced in the spacetime by the gravitational waves as they propagate through the LISA detector are 
measured as the Doppler shifts in the laser frequency. The measured signals will have the various noises along with 
the Doppler shift produced by the gravitational radiation. In the last section we have studied various combinations 
of beams which cancels the laser phase noise and optical bench acceleration noise. In this section we investigate the 
response of the detector for these combinations. We compute the transfer functions for the generators and also their 
linear combinations. The laser phase noise and optical bench acceleration noise is then also cancelled for the linear 
combinations. However, noises such as the shot noise and the acceleration noise of the proof masses do not cancel 
out. In the following subsections we set up the coordinate system adapted to the LISA geometry and then go on to 
compute the response of LISA. 


A. Parameterisation of the interferometer 

Figure [l] describes the LISA configuration. We choose a coordinate system in which the LISA triangle is at rest. 
Although this coordinate system is in motion relative to the usual coordinate systems normally encountered, we will 
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find such a system of coordinates convenient for further analysis, such as computing SNR’s of monochromatic sources 
etc. 

The unit vector w connecting the origin and the source is parameterised by the source angular location ( 9 , </>), so 
that 


w 


/ sin 9 cos 4> 
sin 9 sin <j) 
\ cos 9 


(4.1) 


the transverse plane is spanned by the unit transverse vector 9 and 0, defined by 

~ dw - 1 dw 

9 = 


(4.2) 


89' T sin 0 d(f> 

As the wave propagates through the LISA triangle, the components of the gravitational perturbation can be written 


as 


hij{t, r) = h+(t-w-r) ( 9 t Qj - fafa) + h x (t - w ■ r) ( 9i<f>j + Ojfc ) , 


(4.3) 


where h + and h x are arbitrary functions describing the two GW amplitudes. 

We consider the effect of this perturbation on the light beam traveling between two points A and B. From this we 
obtain the complete response. Let ta and Tb be the position vectors of points A and B respectively. Then the line 
element of the spacetime along the beam null ray obeys, 

0 = dt 2 — dx 2 — dy 2 — dz 2 + hijdx l dx\ (4.4) 

where the i,j run over space indices only. If fi is the unit vector directed from A to B, we have 


dx 1 = n l d \, 


where A is the Euclidean length. The equations (4.4) can be expressed as, 

0 = dt 2 - d\ 2 [1 - hijnV] , 


or equivalently 


From eq. (4.3) we get 


dX = dt 


1 + -h(t 



h{t) = h + (t)£+(8, </>) + h x (f)£ x (9, <t >), 

f £f = {e-nf^^-n ) 2 
\ £x = 2 (^ ' n)(4> ■ n) 


(4.5) 

(4.6) 

(4.7) 

(4.8) 

(4.9) 


B. /i-sensitivity of one arm 

We now apply the above analysis to compute the Doppler response of the laser beam along one arm of the LISA 
detector. Let the beam start at t = to from the point ta and travel towards the point tb and reach it at t = t\. Then, 
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r(t) = r A + (t-t 0 )n, f(ti) = r B . 
The line element along this path satisfies the equation, 


dX = dt 




u> • n)t 


w • VA T t 0 w ■ 



(4.10) 


(4.11) 


The global travel time t\ — to is given by the integral : 


1 f* 1 

L = t 1 — t 0 + — / h[t(l — w ■ n) — w ■ f A + tow ■ h]dt. (4-12) 

2 J to 


It is convenient for many purposes to pursue our analysis in the Fourier domain. We Fourier transform the GW 
amplitude h: 


h(t) = dfl h(Cl)exp(—iflt ), 


(4.13) 


and the travel time can be expressed as, 

t\ — to = L— — J dflh(fl) exp [— if! (1 — w ■ h) t\ exp [iilw • r A \ exp [—ifltow ■ n] dt. 

In the zeroth order of the integral, we have t\ — to = L 1 and we obtain: 

1 


ti — to — L — 


dflh(fl) exp (iflw • f A ) exp (—ittti 


exp (iflLw ■ n) — exp(ittL) 
—*fl(l — w • n) 


(4.14) 

(4.15) 

(4.16) 


The phase change over that time interval is <E> = co(t A — to), where co = 2Ttv opt is the optical circular frequency. We 
can assume that the time t\ is the current time and to the retarded time, so that the phase is <I>(t) = co{t — to): 


<I>(t) = coL — 


co 


2(1 — w ■ n) 
exp {i£lLw ■ n) — exp(iXlL) 
—iQ 


dflh(fl) exp (iflw • r A ) exp (—ittt) 


By taking the time derivative, we get the instantaneous frequency, 

6v{t) 1 ci<b(<) 

Vopt co dt 

In the time domain and using r A + Ln = r B we finally get, 

5v{t) —1 


v opt 2(1 - w ■ h) 

In the Fourier domain we may express this result as: 
5v(£i) h( 0) 


[h (t — w • r B ) — h {t — w • r A — L)\ . 


I^opt 


2(1 — w ■ n) 


exp + w • r A )\ (1 — exp [—if2L(l — w • fi)]) 


(4.17) 

(4.18) 

(4.19) 

(4.20) 

(4.21) 
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C. /i-sensitivity of the elementary Doppler data 


In this section we compute the expression for the transfer function for the six elementary beams given in the eq. 


(2.3). These beams further can be combined with suitable delays as described in the previous sections for achieving 
cancellation of laser phase noise. 

The Doppler shift being expressed in the Fourier domain as, 


Su(n) 

l^opt 


= h+(n)F+(n) + h x (n)F x (n). 


(4.22) 


where F +;X (Q) are transfer functions. We can compute the transfer functions for the combinations Ui, V. We just 
give below F Ul . + x and Fy 1 . + X \ the others are obtained by cyclic permutations: 


F Uy+,x = 


F v i;+,x = - 


e in(w-r 3 +L 2 ) 

2 (1 + w ■ h 2 ) 
e i(l(w-f 2 +L 3 ) 


2 (1 ~w-n 3 ) 


(l - e -ML 2 (l+w ^ 2 + ^ 

(> 


-iQ,L3(l—W-n3 ) 


^3;+, x 5 


where, 


& ; + = (6 ■ hi) 2 - (<f> ■ hi) 2 = 2 (0 ■ fn)((j> ■ hi). 


(4.23) 


(4.24) 


To implement the cancellation of laser phase noise these elementary beams must be combined with suitable time- 
delays. We notice that in the Fourier domain the delay operators get replaced by simple multiplicative factors as the 
following computations show. This is one of the advantages of the Fourier analysis. The delay operators introduced 
in the section are £) such that for any function of time f(t), we have 


Ei*m = at-Li), 


which in the Fourier domain in nothing but 


Ei*f(Sl) = ei (fi)/(fi), 


(4.25) 


(4.26) 


where the are simple factors: 


ei = e 


iQLi 


(4.27) 


Thus operator polynomials in E t become actual polynomials in e* in the Fourier domain. This particularly simple fact 
can be used to advantage for the simple but astrophysically important sources, namely the monochromatic sources 
considered in the next section. 

Below in the figures 2(a) and 2(b) we present the transfer functions Fy 1 for both polarisations: 
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Log(f) 



Log(f) 

FIG. 2. Log Log plot of the psd |FVi. + | 2 and |FVi. x | 2 , of the corresponding transfer function are displayed in (a) and (b) 
repsectively, as a function of frequency for 6 = 0, <j> = 0. 

The other transfer functions show similar characterstics. 


V. MONOCHROMATIC SOURCES 

A. Noise 

We recall that the laser frequency noise and optical bench motion noise can be cancelled by taking appropriate 
combinations of the beams. In this scheme the noises that do not cancel out in the module of syzygies are the 
acceleration noise of the proof masses and the shot noise. These then form the bulk of noise spectrum which we 
obtain below for any given data combination X. We compute the noise power spectral densities for the generators 

XW. 

The beam with the signal and the various noises can be written as, 

U 1 = E 2 C 3 - C{ +2 h 2 -v[ + h v i + Y£\ ot (5.1) 

U 1 = -E 3 C* 2 + Cl + 2n 3 • V, - hv i - YyX ot (5-2) 

Z 1 = {z\ - z{) + h 3 • vi + h 2 • El. (5.3) 
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The other beams can be obtained by taking cyclic permutations. Here V\ and are the random velocities of the 
proof masses, in the left and right branches respectively, in space craft 1 . 

Let the noise cancelling combination X be given by, 

X =p i V i + q i U i + r i Z i , (5.4) 


where the 9-tuple ( pi , qi, ri ) is in the module of syzygies. This combination cancels laser phase noise 
bench acceleration noise, where as the shot noise and proof masses acceleration noise do not cancel, 
we obtain the power spectral density of X for the two noises, 

3 


and the optical 
Using eq. (5.4) 


^ proof mass 


= J2{\2 P i + r l \ 2 + \2q l + r i \ 2 )SP ro °f 


(5.5) 


3 

c x 2 (b,: | 2 +1*| 2 ) s shotnoise , 

i=l 

where S proo ^ mass is obtained from u) and v*. Here following the literature |22j we take, S proo f mass 
10” 48 [//1 Hz]~ 2 Hz- 1 and S shotnoise = 5.3 x 10" 38 [//1 Hz] 2 Hz -1 . 

Explicit simplified expressions for the noise may be obtained by assuming, 

ei = e 2 = e 3 = e lQL 

In the particular cases of the generators X^ A \ A = 1, 2, 3,4 we obtain, 

S xW {f) = [16 sin 2 (2nfL) + 32 sin 4 nfL)] S proofmass + [8 sin 2 (nfL) + 8 sin 2 (27r/L)] S shotnoise 

S xW {/) = 24 sin 2 ( 7 r fL)S proofmass + QS shotnoise 

S xm {f) = [I6sin 2 (7r/L) + 8 sin 2 (37r/L)] S proofmass + QS shotnoise 

S xW (/) = [I6sin 2 (7r/L) + 8sin 2 (37r/L)] S proofmass + QS shotnoise 

The plots of these noise spectra are shown in the figure |[ 


(5.6) 
2.5 x 


(5.7) 

(5.8) 

(5.9) 
(5.10) 



FIG. 3. Log Log plot of the noise spectra for the generators The curves for and 


coincide in the figure. 
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B. Signal 


Monochromatic sources are simplest among the gravitational wave sources. There are a number of important 
objects such as pulsars, rotating neutron stars, coalescing-binaries with sufficiently low mass may be considered as 
emitting monochromatic GW radiation. We will call those sources to be monochromatic which even if they change a 
little in frequency in a given observation time, the fractional change in SNR for the optimal data combination does 
not fall below a preassigned limit. We could take this limit to be few percent, but for concreteness we fix the limit at 
1 %. The observation time T we take to be one year. 

For binary stars, the relevant quantity which decides the evolution in the GW frequency at a given frequency /o is 
the so called chirp mass M. = /r 3 / 5 M 2 / 5 , where /.i is the reduced mass and M is the total mass of the binary system. 
We assume a Newtonian evolution for the binary system which gives the rate of change of GW frequency / as, 

/=“, (5-11) 

O T c 


where r c is the Newtonian coalescence time measured from the epoch when the GW frequency is / 0 and the ‘dot’ 
denotes derivative with respect to time. The r c is given by, 


t c 


2.15 x 10 6 


5 


' M ' 

3 

[ /0 1 



ImHz 


years, 


(5.12) 


where M 0 is the solar mass. 

A limit on the rate of change of frequency can obtained by considering the total change in the frequency A/ during 
the period of observation T. That is, 


A / = fT. 


(5.13) 


Inverting the above equations we obtain a limit on the chirp mass M: 


M < 175 
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r h i 

5 

[ A / 1 

ImHz 


1 [iHz 


M q . 


(5.14) 


In our investigation we take the bandwidth A/ by allowing SNR to change by 1% at the frequency /o- Numerically, 
we estimate A/ for various values of / 0 . The table below shows the upper bound for M. at various frequencies / 0 : 


/o in mHz 

A/ in fi Hz 
for 1% change in SNR 

M 

M @ 

0.1 

1.0 

27705 

1 

9.9 

691 

2 

22 

243 

10 

1130 

74 


Here our goal is to seek a data combination which optimises the SNR for a monochromatic source with given 
polarisation parameters and direction information. A convenient set of polarisation parameters are the angles e and 
if} describing the orientation of the angular momentum vector in space. The direction to the source is described by 
the polar angles 9 and (j> in the coordinate system tied to LISA. 

For a monochromatic source the wave form can be written as 

1 4- cos 2 f 

h+(t) = A --- cos 2 ip cos(2 ut), 

h x (t) = A cose sin2-i/> sin(2wf). 
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(5.15) 

(5.16) 




In the Fourier domain we have. 


1 -f- OOS"^ f 

h+(Q)=A —^- cos 2^, 

h x (fl) = —iA cose sin2 if.’. 


(5.17) 

(5.18) 


The response for the signal at the detector can now be written as, 

3 

hx = ^ [Pi ( Fvi-+h + + Fvi-xhx ) + ft (Fm-,+ h+ + Fui- X h x )\ , (5.19) 

i= 1 

where p^s and qi s are in the module of syzygies. From eq. 
for the generators and it can we written as, 

3 

S x (f) = J2 [(1 2 Pi + ^\ 2 + 1 2ft + n\ 2 ) S proofmass + (| K | 2 + |ft| 2 ) S shot ”°“ e ] . (5.20) 

i=1 

Expression for the signal to noise ratio (SNR) for the monochromatic source simplifies to, 


(5.5) and (|5.f|) we can compute the total noise spectrum 


SNR = 


\hx\ 1 
y/SxU) J ' 


(5.21) 


We plot the sensitivities of the generators X as function of / in figures |^(a) and||(a) by fixing the angles 6 and </>. 
It is also important to understand the angular dependence of the sensitivity of generators, which are plotted in the 
figures ||(a)-||(d) at a frequency of / = 1 mHz. The sensitivity S is defined following |2^| , 


S = 


\hx\ 


(5.22) 


where B = 1/T, where T is the observation time which we take to be one year. The number 5 corresponds to an SNR 
of 5. 
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FIG. 4. Plots of log S of the generators X^ A \ A = 1, 2, 3, 4, are displayed in (a), (b), (c) and (d), respectively, as function 
of 6 and tj> for f = 1 mHz and SNR = 5. 


C. Maximisation 


In this subsection our goal is to maximise the SNR for a given monochromatic source over the set of noise cancelling 


combinations. These combinations can be generated by the generators given in the eq. ([T0|) and (3.7). The SNR 
corresponding to the each of the generators (X^ A \ A = 1 to 4) as a function of frequency are shown in the figure. 
However one must maximise the SNR over an arbitrary linear combination of X^ A \ This goal is difficult to achieve 
since it involves a maximisation over a space of six tuples of polynomials which is essentially a function space. In 
order to make the problem tractable and still achieve adequate results we restrict the polynomials to be constants. 
This approach does not fully optimise the SNR but it comes quite close to the optimal solution. Our approach can 
be thought of as a zero’th order approximation. 

A linear combination of the generators can be written as 


x = J2 a (A) X(A \ 


(5.23) 


A— 1 


here, ol^a) (f° r A= 1 to 4) are a set of real numbers. Since a scalar multiple of X will not yield anything new, we set 
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one of the cn’s, say, cqi) = 1. Thus the SNR now becomes a function of three parameters ayyi = 2,3,4, which are 
just real numbers and our objective is to maximise the SNR with respect to a^y 

In order to carry out the analysis efficiently and elegantly we find that it is useful to define complex noise vectors 
N pertaining to X b 4 ) as follows: 


«"•> = [ys[(2 V ^+r\ A \yw 1 (i q f' +r^),ys;^\s/w 2 ,^) 


(5 24) 


where, p\ A \ qA and r ? -' 4i corresponding to generators X^ are given in the eq. (3.6) and ( 3.13 ) and 

Si = SP r °°f mass anc i s 2 = cjshotnoise _ have N^ £ C 12 the 12 dimensional complex space and the 

usual scalar product C 12 induces a norm; ■ N^* = ||7V ^' 4 ^|| 2 gives the noise psds corresponding to the basis 

In a similar fashion one can also write the signal corresponding to a particular basis element. We first define the 
polynomial 6 -tuple for each generator X ^ as follows: 


p (A) = (pAA A) ) , 

and the GW signal 6 -tuple as, 

H* = ( F Vi -+h + + F Vi - tX h x , Pu,;+h+ + F Ui;x h x ). 

The signal for a specific generator X ( - A> is then written as, 

h (A) = P (A) .H* 

and the corresponding SNR is given by, 


SNR( a) = 


IM A )I 


(p<w • p*)(p( A ) • H*y 


(5.25) 


(5.26) 


(5.27) 


(5.28) 


||7Vh4)|| V NW-NW* 

For an arbitrary linear combination X (eq. ( 5.23| )) the noise vector and the signal vector can be expressed as 

N = a (A) N( A \ P = a {A) pW, (5.29) 

where summation convention has been used. The signal is just the scalar product h = P.H* = We omit 

subscripts X on these quantities. 

In this notation the SNR of the combination ([5.2S ) can be written as, 


SNR iiJir 

Writing out explicitly the sums in the scalar products, 

c*(A )a(B)h (A) h^ 


Csnr y = 


a (A) a (B) N( A ) -NW* 


(5.30) 


(5.31) 


Maximisation with respect to cc( 2 ), 0 ( 3 ), 0 ( 4 ) leads to the following three conditions which must be obeyed by qw in 
order to yield the maximum SNR for X: 


5ft(/i/i«*) 3ft(AW«*) 


\h \ 2 


2 ’ 


(5.32) 


where 3?(x) denotes the real part of the quantity x. 
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To demonstrate the usefulness of the formalism, we consider just two generators X^ and . We take the 
«(!) = 1, «( 2 ) = a and other two a s zero. Then, 


X = X^+aX^K 


The eq. (5.31) reduces to the form, 


where, 


a\ + 2b\a + c\a 2 
a 2 + 2 b 2 a + c 2 a 2 ’ 


(5.33) 


(5.34) 


a 1 = \h 1 \ 2 , b 1 = M(h 1 h 2 *), ci = \h 2 \ 2 , 
a 2 = \N 1 \ 2 , b 2 = $t{N 1 N 2 *), c 2 = \N 2 \ 2 . 


(5.35) 


The condition for the optimisation (5.32) simplifies to 

(6ia 2 — ai6 2 ) + (cia 2 — aic 2 )a + (ci& 2 - 6iC 2 )a 2 = 0. (5.36) 

The two roots of the eq. (|5.3G|) can be obtained. Here, a is a function of the parameters /, 9, </>, e and ip. One of 


the solutions of the eq. (5.36) correspond to the maximum and other correspond to the minimum of the SNR. In a 


similar fashion one can maximise the SNR by taking any two of the four generators given in the eq. (3.6) and by 
taking appropriate avp in the eq. ( |5.23 ). We have seen in several cases that maximising over just two generators 
yields remarkably good results. 



FIG. 5(a). Log Log plot of the sensitivity S for the generators X ^ as function of / for 9 = 0, (j> = 0 over one year observation 
period for SNR = 5. The curve [X1,X2,X3,X4] depicts the sensitivity of the linear combination of four generators which 

gives maximum SNR. 
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FIG. ||(b). Plot of coefficients dm which gives the maximum SNR for linear combinations of all the four X^ A ' 1 as 
function of / for 9 = 0 and (f> = 0. 



FIG. 6(a). Log Log plot of the sensitivity S for the generators X^ as function of / for 6 = 7t/4, (j) = 7r/4 over one year 
observation period for SNR = 5. The curve [X1,X2,X3,X4] depicts the sensitivity of the linear combination of four generators 
which gives maximum SNR. 
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FIG. |](b). Plot of coefficients a(A) which gives the maximum SNR for linear combinations of all the four X^ 4 ) as 
function of / for 9 = 7r/4 and <j> = 7t/4. 


This simple case demonstrates that one can use the solutions given in the eq. (3.6) to get a better SNR. However, 
to get full advantage one needs to maximise the SNR over the three a’s. In order to optimise the SNR given by the 
general combination we resort to numerical methods since there is no straight forward method for solving the coupled 
algebraic equations given by ( 5.32 ). We use the Powell’s method as given in Jg| for maximising the SNR over the 
parameters atA)- The sensitivity S for the generators X^ A> and for the maximal SNR combination of X^’ s denoted 
by [X1,X2,X3,X4] as a function of frequency / has been plotted in figures |](a) and |](a). The coresponding values of 
oi(A) are shown in figures |b) and ib). 


VI. CONCLUDING REMARKS 

We have presented in this paper a rigorous and systematic procedure for obtaining data combinations which cancel 
the laser frequency noise based on algebraic geometrical techniques and commutative algebra. The data combinations 
cancelling the laser noise have the structure of a module called the module of syzygies. The module is over a ring of 
polynomials in three variables, corresponding to the three time delays along the three arms of the interferometer. Our 
formalism can be extended in a straightforward way to include (cancel) the Doppler shifts due to the motion of the 
optical benches. This module provides us with a choice of data combinations which are in turn linear combinations 
of the generators of the module. We use this parametrisation to maximise the SNR over frequency for one class of 
GW sources, namely, those that are monochromatic. We observe that in the plot of sensitivity verses direction angles 
for the generators X^ A \ namely figures [| (a)-(d), the sensitivity has several peaks. It may be possible to employ 
this property to optimally resolve binaries in the confusion noise regime by considering suitable data combinations 
which would be sensitive to specific directions in the sky. Also here we have investigated monochromatic GW sources. 
We believe that this formalism may also be applied successfully to other type of GW sources eg. stochastic GW 
background. 
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APPENDIX A: GENERATORS OF THE MODULE OF SYZYGIES 


We require the 4-tuple solutions (p 3 , 51 , 52 , 93 ) to the equation: 

(xyz - l)p 3 + (xz - y)qi + x(l - z 2 )q 2 + (1 - x 2 )q 3 = 0 , 


(Al) 


where for convenience we have substituted x = E\,y = E 2 , z = £ 3 . p 3 , gi, q 2l q 3 are polynomials in a;, y , 2 with integral 
coefficients i.e. in Z[x,y,z]. 

We now follow the procedure in book by Becker et al. [Q . 

Consider the ideal in Z[x , y , z ] (or Q[x, y, z ] where Q denotes the field of rational numbers), formed by taking linear 
combinations of the coefficients in eq.(Al) f\ = xyz — 1, f 2 = xz — y, f 3 = x(l — z 2 ), fi = l — x 2 . A Groebner basis 
for this ideal is: 


G = {gi = z 2 - 1, 92 = y 2 ~ 1,23 = x - yz} . (A2) 

The above Groebner basis is obtained using the function GroebnerBasis in Mathematica. One can check that both 
the fi,i = 1,2,3,4 and gj,j = 1,2,3 generate the same ideal because we can express one generating set in terms of 
the other and vice-versa: 


fi — dijgj , 


9 3 ~ c jifi 


where d and c are 4x3 and 3x4 polynomial matrices respectively, 


and are given by, 


^ 1 z 2 yz ^ 

y 0 2 

—x 0 0 

V -1 - 2 2 -(x + yz) ) 



0 

-y 

Z 


-x z 2 — 1 \ 

0 0 1 
1 0 ) 


(A3) 


(A4) 


The generators of the 4-tuple module are given by the set A (J B* where A and B* are the sets described below: 

A is the set of row vectors of the matrix I — d.c where the dot denotes the matrix product and I is the identity 
matrix, 4 x 4 in our case. Thus, 


ai = (1 - 2 2 , 0 , x — yz, 1 - 2 2 ) , 
a 2 = (0, 2(1 - 2 2 ), xy - 2, y(l - z 2 )), 
a 3 = (0, 0,1 - x 2 ,x(z 2 - 1)) , 

04 = (z 2 , xz, yz, z 2 ). (A5) 


We thus first get 4 generators. The additional generators are obtained by computing the S-polynomials of the Groebner 
basis Q. The S-polynomial of two polynomials 51,52 is obtained by multiplying 51 and 52 by suitable terms and then 
adding, so that the highest terms cancel. For example in our case 51 = z 2 — 1 and 52 = y 2 — 1 and the highest terms 
are z 2 for gi and y 2 for g 2 . Multiply 51 by y 2 and g 2 by z 2 and subtract. Thus, the S-polynomial p 32 of 51 and g 2 is: 

P 12 = y 2 gi - z 2 g 2 = z 2 -y 2 . (A6) 

Note that order is defined (x » y » z) and the y 2 z 2 term cancels. For the Groebner basis of 3 elements we get 
3 S-polynomials Pi 2 ,Pi3,P23- The Pi :j must now be rexpressed in terms of the Groebner basis Q. This gives a 3 x 3 
matrix b. The final step is to transform to 4 -tuples by multiplying b by the matrix c to obtain b* = b.c. The row 
vectors b*, i = 1 , 2 , 3 of b* form the set B *: 
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bt = ( 1 - * 2 , y(z 2 - 1), *(1 - y 2 ), (y 2 - l)(z 2 - 1)), 
b * 2 = (0, z( 1 - z 2 ), 1 - z 2 - x(x - yz ), (a: - yz){z 2 - 1)) 
b * 3 = {x- yz, z-xy, 1 - y 2 ,0). 


(A 7 ) 


Thus we obtain 3 more generators which gives us a total of 7 generators of the required module of syzygies. 


APPENDIX B: MATRICES OF CONVERSION BETWEEN GENERATING SETS 

In this appendix we list the three sets of generators and relations among them. We first list below a, ( 3 , 7, 

a = (1 ,z,xz, 1 ,xy,y), 

( 3 = (xy,l,x,z,l,yz), 

7 = {y,yz, l,xz,x,l), 

C = {x, y, Z, x, y,z ). (Bl) 

We now express the a,; and b * in terms of a, ( 3 , 7, 

ai='i-zC, 
a 2 = a — z/ 3 , 

a 3 = — za + f 3 — x'y + xzC,, 

(I4 = zC, 

bl = -ya + yz (3 + 7 - zC, , 

£>2 = (1 — z 2 )/3 — £7 + > 

b * 3 = ( 3 -y<;. (B 2 ) 

Further we also list below a, ( 3 , 7, £ in terms of X^: 

a = XW , 

13 = I (4) , 

7 = -X (1) + zX (2) 3 , 

c = X (2) . (B 3 ) 

This proves that since the a*, b* generate the required module, the a, ( 3 , 7, ( and X^ A \A = 1 , 2 , 3,4 also generate the 
same module. 

The Groebner basis is given in terms of the above generators as follows: G ^ = £, G^ = X^\ G^ = (3 , G = a 
and G^ = 03. 
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